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Abstract 

We evaluate the cross section of the process e + e~ — ► 7^(3872) in terms of the 
content of the D*°D° 'molecular' component in the wave function of the resonance 
X(3872). If this component is dominating, the cross section of the reaction e + e~ — ► 
7^(3872) can reach up to about 10 -3 of that for e + e~ — ► D*D* at energy slightly 
above the D*D* threshold, and the considered process can be a realistic source of the 
X(3872) particles for the studies of this resonance. 



The extreme proximity of the mass of the X(3872) resonance fl] to the D°D*° threshold, as 
well as the co-existence of the decays X(3872) — > 7i + 7i^J/ip and X(3872) — > n + n~ n° J / 2J, 
strongly suggestsj3IUI3ini a significant presence of a 'molecular' [7 a D°D*° + D*°D° compo- 
nent in the wave function of X(3872). Such unusual structure of this hadronic state makes 
further studies of its properties very promising for understanding the strong dynamics of 
heavy mesons. At present the X(3872) resonance is observed experimentally only in the 
decays of B mesons B — > X K [HI Ej and in inclusive production in proton - antiproton col- 
lisions at the Tevatron|10[ ITT] . Both these types of processes are quite rare and also present 
significant challenges for precision measurements of the parameters of the discussed reso- 
nance. In particular^, neither the total width of X(3872) is yet resolved (the current limit 
is Tx < 2.3 MeV), nor its mass is known with a precision sufficient to determine the binding 
energy of the molecular component (the current average value Mx = 3871.2 ± 0.5 MeV coin- 
cides within the errors with M(D°) +M(D*°) = 3871.2 ±0.8 MeV). In this paper we consider 
the process e + e" — » 7X(3872) at the cm. energy within few MeV of the D*°D*° threshold, 
where the kinematical simplicity of the process would hopefully allow more detailed studies 
of X(3872). We estimate that the cross section a[e + e~ — > 7X(3872)] is likely to be at least 
about 10 -3 of the cross section for the production of D*D* meson pairs, i.e. in the range 
of about 1 pb, which makes realistic a study of the discussed here process in a dedicated 
experiment. Moreover, the energy dependence of the cross section is sensitive to the binding 
energy of the molecular component. Thus a study of this dependence can provide a better 
accuracy of determining the mass of X(3872) relative to the D°D*° threshold than a direct 
mass measurement. 

In order to estimate the cross section of the discussed process we calculate the absorptive 
part of the production amplitude due to the process e + e~ — > D*°D*° — ► 7^(3872) with 
on-shell D* mesons, as given by the unitarity relation. We find that this contribution to 
the amplitude is rapidly changing with the cm. energy of the process. The contribution of 
other intermediate states to the amplitude, which potentially could destructively interfere 
with the calculated amplitude is a slowly varying function of energy, so that such destructive 
interference cannot occur at all energies in the considered range. Thus the value of the cross 
section can be estimated as being at least that given by the calculated part of the amplitude. 

Proceeding to the calculation of the absorptive part of the amplitude we write the Fock 
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decomposition of the wave function of X(3872) in the form 

D D*° + D*°D „ , 

*x = a 7= h l^diipi , (1) 

V2 i 

where the molecular component is explicitly separated, and the sum runs over 'other' states. 
The crucial difference between the molecular and 'other' components of X(3872) is that 
due to a very small binding energy w the neutral D and D* mesons mostly move at long 
distances, beyond the range of the strong interaction, while the 'other' states are localized 
at shorter distances typical of the strong interaction. At long distances the coordinate wave 
function of the meson pair is that of a free S-w&ve motion and is proportional to exp(— nr)/r, 



where k = \J1 m r w ~ 44 MeV J w(MeV) is determined by the reduced mass m r « 966 MeV 
in the D°D*° system and the binding energy w. With the coordinate wave function of 
this component normalized to one, the coefficient a$ in the expansion (P) determines the 
statistical weie ht |a | 2 of the molecular state (D°D*° + D*°D°)/y/2 in the wave function of 
X(3872). The notion of the resonance X being mostly a molecular system corresponds to 
this statistical weight factor of order one. The numerical value of | a, 1 2 is presently unknown, 
in a model calculation ^2] this weight factor is estimated as 0.7 - 0.8 at w — 1 MeV. 

In what follows we calculate the contribution of the 'peripheral' D°D*° + D*°D° com- 
ponent of the X(3872) resonance to the absorptive part of the amplitude of the process 
e + e~ — > D*°D*° — > 7X(3872), where the latter transition proceeds due to the underlying 
radiative decay D*° — > D 0/ y (D*° — > D 0r y), in analogy with the previously discussed [3 EI] 
decays of the X resonance, X DDj. The amplitude of the radiative decay of the vector 
meson has the form 

A(D*° — > J D°7) = /x e ijk e? kj a k , (2) 

with a and e being the polarization amplitudes of the vector meson and the photon, and 
k being the photon momentum. 1 The vector meson decay rate Tq = T(D* — > D^) is then 
related to the transition magnetic parameter /x as 



3tt 



(3) 



where Uq ~ 137 MeV is the photon energy in the decay. Numerically the decay rate can be 
estimated from the datapP as T = 26 ± 6KeV. The amplitude of the transition D*°£)* — > 
lr The non-relativistic normalization for the heavy meson states is used throughout this paper. 
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7^(3872) due to the 'peripheral' component of the X can then be written in terms of the 
momentum-space wave function (f)(q) of this component: 



A(D*°D*° ^ 7 I) 









a k ( 







P+ 2 



(4) 



where, in addition to the notation conventions in Eq.(J2J), b stands for the polarization am- 
plitude of the initial D* meson and x is the polarization amplitude of the produced X 
resonance. The relative minus sign between the two terms originating from the amplitudes 
of the processes D*° — > D 0/ y and D*° — > D°7 is due to the opposite C parity of the X 
(J p c = 1++) and of the photon HHJ- 

In the present calculation we describe the peripheral component by a wave function of a 
free-motion with an ultraviolet regularization at large momenta |14t I13j : 



4>(q) = V87TKC 



q 2 + k 2 q 2 + A 2 
where the normalization constant c is given by 



(5) 



A (A + k) 

A- K 



(6) 



and the regularization parameter A is determined by the inverse size of the strong interaction 
region in the X resonance. The obvious reason for introducing the cutoff A is that the free- 
motion description of the meson pair inside X is applicable only at distances beyond the 
range of the strong interaction and such description generally fails at short distances, where 
the mesons overlap and the meson pair strongly mixes with the 'other' states in the expansion 
(HJ). Thus introducing the parameter A is a way of explicitly separating the 'peripheral' part 
from the 'core'. 

The other ingredient in the calculation of the amplitude of the process e + e _ — > D*°£)* — > 
7AT(3872) is the production amplitude for the D*D* pair in e + e~ annihilation. At energy 
E = 2M(D*°) + W near the threshold this amplitude can be generically written in the form 



A{e + e 



*o 7-1*0^ 



D *u D 



A (l-p)(S-b)* + 



2^5 



MjiPk 



di b k + a k bi - - 5 ik (a ■ b) 



(7) 



where j = (eje) stands for the current of the incoming electron and positron, p is the 
momentum of one of the mesons (D*° for definiteness) in the cm. frame, and Aq and A2 
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are the factors corresponding to production of the vector meson pair in the states with 
respectively the total spin S = and S = 2. It can be also noted that the amplitude in 
Eq.(|7|) describes the production of mesons in the P wave. Another kinematically possible 
amplitude, the F-wave, should be small near the threshold, i.e. at a small W. Both A and 
A 2 are generally functions of the excitation energy W. Furthermore, their dependence on the 
energy near the threshold is known to be nontrivial due to the ■0(4040) resonancepQ, with 
possible further complications in the immediate vicinity of the threshold fTHj . Neither 
the relative magnitude nor the relative phase of the amplitudes Aq and A 2 is presently 
known, but both of these can be measured from angular correlations [17^ These amplitudes 
determine the total cross section for production of D*°D*° in e + e~ annihilation: 

a(e + e- ^D*°D*°)= [ \A(e + e~ — > D*°D*°)\ 2 2ty 5 (w — —\ ^P =C H^. (\A \ 2 + \A 12 
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where m = M(D*°), p = \p\, and C is an overall constant related to the average value of 
the current \j\ 2 . The specific value of the latter constant will not be essential in further 
calculation, since it cancels in the ratio of the cross sections. It should be pointed out that 
the nonrelativistic expression for the phase space is used in Eq.fjHJ) corresponding to the 
nonrelativistic normalization of the states of the heavy mesons. 

The discussed here absorptive part of the amplitude of the process e + e~ — > 7X(3872) 
due to the D*°D*° intermediate state is found from the unitarity relation and the amplitudes 
(J1J and (0) in the standard way: 

^Abs(e + e- -> jX) = (9) 
\JY, A ( e+e ~ D*°D*°) A(D*°D*° -> 7 X) 2tt 5 ( W - ^ ^ P 



pol 

fiaopm [~k 
2a; 2 V n 



F ^ijk kj 



m) (2tt) 3 

xlO-k) (A -^\+j k (k-r)^=A 2 



In the latter expression u — \k\ is the energy of the photon, the sum goes over the polar- 
izations of the vector mesons in the intermediate state, and F stands for the dimensionless 
form factor: 

l r 1 _ a , / ' k 



J (P-k)4> [p-2) dcos6 ' (10) 



y/SlTK J-l 

where 9 is the angle between the vectors p and k. Using the expression (J3J) for the 'peripheral' 
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wave function, one readily finds the form factor as 



F=^- 

puj 




(p-cj/2) 2 + /t 2 ^ 4 / (p-cj/2) 2 + A 2 

with the normalization coefficient c given by Eq.(|5jl. 

The absorptive part of the amplitude in Eq.Q corresponds to the cross section 

d 3 k 



(11) 



<x A bs(e + e- _> D*°f)*° ^ 1 X) = J |A Abs (e+ e - -> 7 X)| 2 27r<5(^ - |fc|) 
'a*| 2 I°o| 2 P 2 m 2 uj kF 2 



C 



12tt 2 




(27t) 3 2cj 

(12) 



where the overall constant C in the latter expression is the same as in Eq.fjHJ). Thus using 
also Eq.© one finds the formula for the ratio of the cross sections: 

<XAbs(e+e- - D*°D*° -> 7 X) 2 r mcu«; 2 |A - A 2 /V5| 2 + (9/20) \A 2 \ 2 

~\ a °\ o. .3 ~ * 1 /t 19 i I /I 19 • (l 3 ) 



a(e+e- -> D*°D*°) 1 U| 2u;gp 

The so-defined cross section o"Abs is (most likely) not the actual value of the cross section, since 
the amplitude of the process e + e~ — > 7X can receive contribution from other mechanisms. 
Nevertheless it is instructive to examine the numerical value and the behavior with energy of 
this quantity as given by Eq.(fT3*]). The dependence on the cm. energy of the factor (k/p) F 2 
is shown in Fig.l for two representative values of the 'molecular' binding energy in X(3872), 
w = 1 MeV (k w 44MeV) and w = 0.3 MeV (k ps 24MeV). This factor peaks at the energy 
where p ~ uj/2 ps 70 MeV. The appearance of this peak is easily understood qualitatively: 
at p ps k/2 the D° meson emerging from the emission of the photon in D*° — > D Qr j moves 
slowly relative to the D*° and forms a loosely bound state. 2 The width of the peak is clearly 
determined by the parameter k. 

As is seen from the plots of Fig.l the numerical value of the factor (n/p) F 2 near its peak 
is of order one. Another factor in Eq.()13|). Vomu/fiuQ) ~ rom/(2c<j 2 ) ps 1.5 x 10~ 3 , sets the 
overall scale of the discussed cross section. The factor in Eq. fjrHj) depending on the presently 
unknown ratio of the (generally complex) amplitudes A /A 2 , takes values between 0.34 (at 
A /A 2 ~ 0.68) and 1.31 (at A Q /A 2 ~ 1-47), and can thus be considered as being of order 
one. Finally, the statistical weight factor |ao| 2 , as discussed, is likely to be large fraction of 
one. Summarizing these numerical estimates, the value of the ratio in Eq. (jl3|) at the peak 
can be estimated as being of order 10~ 3 , although the uncertainty is presently large. 
2 The same situation arises at p s» — k/2 for the D° meson emerging from D*° — ► D ^. 
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Figure 1: The factor nF 2 /p vs. the excitation energy W above the D*°D*° threshold at 
representative values of the binding energy w in X(3872) and the ultraviolet cutoff parameter 
A: w = lMeV, A = 200 MeV (solid), w = 1 MeV, A = 300 MeV (dashed), w = 0.3 MeV, 
A = 200 MeV (dashdot), and w = 0.3 MeV, A = 300 MeV (dotted). 

In absolute terms, the measured |1 5 j cross section a(e + e~ — ► £)*°|)* ) at E — 4015 MeV, 
i.e. at the energy above the D*°D*° threshold W ~ 1.6 MeV is about 0.15 nb. This cross 
section grows from the threshold as p 3 . With this factor taken into account the peak of the 
quantity o"Abs(e + e~ — ► jX) shifts to a slightly higher value of p, p ~ 100 MeV, correspond- 
ing to W Rj 5 MeV, where according to Eq. (|T2"|) and the presented estimates, it should be 
numerically of the order of 1 pb. 

The considered mechanism of the process e + e~ — > describes a 'soft' production of 
its peripheral D°D*° + D*°D° component in radiative transitions from slow £)*°]J* pairs. 
Generally, one can also expect a presence of states with charged mesons, D + D*~ + D~D* + , 
within the X resonance. Therefore a contribution of the mechanism e + e~ — > D* + D*~ — > ^X 
to the considered here process merits discussion. However, this contribution in fact should 
be very small for at least three reasons The mass of the pair of charged mesons is 
by 5 ~ 8 MeV heavier than that of the neutral ones. For this reason the charged mesons 
in this component are separated by shorter distances: the corresponding parameter k is 
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approximately y/2m r 5 ~ 125 MeV, and as a result |T3| the statistical weight of such state 
should be significantly smaller than for the pair of neutral mesons. Furthermore, the C- 
conjugate processes D* + — > D +r y and D*~ — > D _ 7 destructively interfere in D* + D*~ — > 
7X (cf. the minus sign between the two terms in Eq.(J3])). The (negative) interference is 
enhanced for the more closely separated charged mesons in X(3872). Finally, the transition 
magnetic coupling fi is noticeably weaker for the charged mesons than for the neutral ones: 
T(D*+ -> D+7) = 1.5 ± 0.5 KeV. 

Other intermediate states with charmed meson pairs, i.e. DD and DD* (DD*), can 
potentially contribute to the discussed process e + e~ — > jX. Indeed, the suitable final state 
arises in the chain e + e~ — > DD — > 7 (DD*+D*D) through the radiative transition D — > 7D* 
( D -> 7-D*) as well as in the chain e+e~ -> DD* + DD* -> 7 (DD* + D*D) through an 
elastic emission of a photon by D* (D*). However, one can readily see that in either of 
these processes the charmed meson emerging after the emission of the photon is very far 
off the mass shell in the scale of k. Thus neither of these processes can proceed due to the 
long-distance peripheral component of the X(3872) resonance, but rather is determined by 
the short- distance dynamics of the 'core' of X. For this reason these contributions, as well as 
other possible mechanisms related to the 'core' dynamics, should be smooth functions of the 
cm. energy on the scale of few MeV around the D*°D*° threshold, where the amplitude given 
by Eq. ([9|) experiences a significant variation. Therefore even under the most conservative 
(and quite unlikely) assumption that these mechanisms cancel the contribution of the latter 
amplitude near its maximum, such cancellation cannot take place at all energies in the 
considered range. Thus the cross section of the process e + e~ — > 7X at an energy within few 
MeV of the D*°D*° threshold has to be at least as large as the above estimates for o"Abs near 
its maximum i.e. of the order of lpb. The latter is a conservative estimate, since we cannot 
exclude that the contribution of those 'other' mechanisms exceeds the calculated amplitude 
and that the actual cross section is larger than dAbs- 
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